Casimir interaction of dielectric gratings. 
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We derive an exact solution for the Casimir force between two arbitrary periodic dielectric gratings 
and illustrate our method by applying it to two nanostructured silicon gratings. We also reproduce 
the Casimir force gradient measured recently [l| between a silicon grating and a gold sphere taking 
into account the material dependence of the force. We find good agreement between our theoretical 
results and the measured values both in absolute force values and the ratios between the exact force 
and PFA predictions. 



INTRODUCTION 

The availability of experimental set-ups that allow ac- 
curate measurements of surface forces between macro- 
scopic objects at submicron separations has recently 
stimulated a renewed interest in the Casimir effect. In 
1948 H. B. G. Casimir showed that two electrically neu- 
tral, perfectly conducting plates, placed parallel in vac- 
uum, modify the vacuum energy density with respect to 
the unperturbed vacuum|2j. The vacuum energy density 
varies with the separation between the mirrors and leads 
to the Casimir force, which scales with the inverse of the 
forth power of the mirrors separation L. 

The Casimir force is highly versatile and tailoring it 
could potentially be useful in the design and control of 
micro- and nanomachincs. While the material depen- 
dence of the Casimir force has been thoroughly studied 
between two plane mirrors (see e.g. 0, US @]), for 
most other geometries exact calculations exist only for 
perfectly reflecting boundaries (see e.g. (?]])• If material 
properties are taken into account, the shape dependence 
of the Casimir force is usually treated using the proximity 
force approximation (PFA) which amounts to summing 
up contributions at different distances as if they were in- 
dependent. 

In a recent paper [l|, Chan et al. present the first mea- 
surement of the Casimir force between a silicon grating of 
high aspect ratio and a gold sphere and demonstrate the 
violation of PFA in this geometry. Corresponding calcu- 
lations taking into account the periodic structure beyond 
PFA, but only for perfect mirrorsjl], turn out to lead to 
a too large deviation from PFA[l|. 

In this Letter we present the first exact calculation of 
the Casimir force between gratings of arbitrary periodic 
structure, where we take explicitly into account the (ar- 
bitrary) dielectric permittivity of the material. We first 
present formulations for the Casimir energy between two 
periodic dielectric gratings and outline the derivation of 
these formulae. We then apply our formulation to the 
situation of two rectangular silicon gratings and show 
that our calculation yields deviations of the real force 
from the PFA prediction up to 24 percents. We also per- 
formed calculations corresponding to the measurement 



by Chan et al. allowing therefore a first quantitative 
theory-experiment comparison. The result taking into 
account the finite conductivity gives a smaller deviation 
of the exact force from the PFA prediction than the cal- 
culation for perfect mirrors. 



GENERAL PROCEDURE 

We consider two periodic dielectric gratings of arbi- 
trary form separated by a vacuum slit. The special case 
of lamellar (or rectangular) gratings is depicted in FigQ] 
The geometrical parameters are the corrugation depth a, 
the period d and the gap d\ . The gaps of both gratings 
are separated by a distance L. For simplicity we will sup- 
pose the space between the two gratings to be filled with 
vacuum with e = fx = 1. 

The physical problem is time and z invariant, so elec- 
tric and magnetic fields can be written in the form: 

Ei(x,y,z,t) = Ei(x,y) exp(ik z z - iut), (1) 
Hi(x,y,z,t) = Hi(x,y)exp(ik z z-iu>t). (2) 

Let us first suppose the upper grating to be absent and 
consider a generalized conical diffraction problem on the 
lower grating. The longitudinal components of the elec- 
tromagnetic field outside the corrugated region (y > a) 
may be written by making use of a generalization of the 
Rayleigh expansion for an incident monochromatic wave: 
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with an integer p. The sums are performed over all in- 
tegers n. All other field components can be expressed 
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FIG. 1: Rectangular gratings geometry. 



via the longitudinal components E z , H z . This solution is 
valid outside any periodic one-dimensional structure. 

We now have to determine the coefficients Rnp , Rn P ' 
for a specific periodic geometry profile. For this purpose 
we rewrite the Maxwell equations inside the corrugation 
region < y < a in the form of first order differential 



equations, ^ = MA, where M is a square matrix of 
dimension 87V + 4, A T = (E z , E x , H z , H x ) and 2N + 1 is 
the number of Rayleigh coefficients considered in every 
Rayleigh expansion. For a rectangular dielectric grating 
the matrix M is a constant matrix. At y — the solution 
has to satisfy the following expansions, valid for y < 0: 
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We then determine the unknown Rayleigh coefficients by 
matching the solution of equations ^ = MA inside the 
corrugation region with Rayleigh expansions ©j© a t 
y = a and expansions ([7 ]) .([8" ]1 at y = 0. Everywhere in 
the calculations we assumed fj, = 1. 

The fields E z and H z are not decoupled for k z ^ 0. 
This is why the reflection matrix Ri for a reflection from 
a lower grating can be defined as follows: 
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Performing a change of variables y = —y' + L, x = 
x' — s (s < d) in (JSj), (J3J, it is possible to obtain the re- 
flection matrix i?2«p for the reflection of an upward wave 
from a grating with the same profile turned upside-down, 
displaced from the lower grating by Aa; = s, Ay = L. 
Note that for the upper grating in FigQ]the special case 
s = is depicted. 

Up to now we considered a diffraction problem on a sin- 
gle grating. In Q the Casimir energy between two bodies, 
the diffraction properties of which can be described by a 
scattering matrix, has been derived in plane geometries 
on the basis of canonical quantization. Roughness correc- 
tions were derived on the basis of a scattering approach 
in [Io| . The path integral method was used to obtain 
multipole expansion of the Casimir energy between the 
two compact objects llj, exact results in spherical ge- 
ometries [ll|, E2] were also derived. 

We outline a novel derivation here, which can be ap- 
plied to various Casimir systems. To obtain the Casimir 
energy we need to determine the eigenfrequencies of all 
stationary solutions of the generalized diffraction prob- 
lem of subsequent diffraction of the electromagnetic field 
on two periodic gratings separated by a gap-gap distance 



L. These eigenfrequencies can be summed up by making 
use of an argument principle, which states: 

1 f d 

-j— j> ln/(w)tL; = ^(Hw )-X^( w °°)' ( n ) 

where luq are zeroes and ujoc are poles of the function 
f(ui) inside the contour of integration. Degenerate eigen- 
values are summed over according to their multiplicities. 
For the Casimir energy we have <f>(uj) = ftw/2. The equa- 
tion for eigenfrequencies of the corresponding problem of 
classical electrodynamics is f(ui) = 0. 

Consider first the plane-plane geometry when two di- 
electric parallel slabs (slab 1: y < 0, slab 2: y > L) are 
separated by a vacuum slit (0 < y < L). In this case 
TE and TM modes do not couple. The equation for TE 
eigenfrequencies is: 



f(oj) = 1 - r 1TE {oj)r 2 TEu P {^) = 0. 



(12) 



Here r\TE{u) is the reflection coefficient of a downward 
plane wave which reflects on a dielectric surface of slab 1 
at y = 0, while r2TEup{^>) is the reflection coefficient of an 
upward plane wave which reflects on a dielectric surface 
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of slab 2 at y = L. One can deduce from Maxwell equa- 
tions that r 2 TEu P (u) = t 2 te{ i -o) exp(2ik y L) (r 2 TeM is a 
reflection coefficient of a downward TE plane wave which 
reflects on a dielectric slab 2 now located at y < 0) . From 
(fT2"| and the analogous equation for TM modes one im- 
mediately obtains the Lifshitz formula by making use of 
the argument principle (jlip . 

For two periodic dielectrics separated by a vacuum slit 
one has to consider a reflection of downward and upward 
waves from a unit cell < k x < 2w/d. Due to the struc- 
ture of the surface, TE and TM modes do not decouple 
anymore, but they are coupled by the diffraction process. 
The equation for normal modes states: 

Ri(u>i)R 2up (uJi)ipi = tjii, (13) 

where tpi is an eigenvector describing the normal mode 
with a frequency Wj. Instead of equation (|12| one obtains 
from (|13p the following condition for eigenfrequencies: 

det(I - i?iHi?2u P M) = 0- ( 14 ) 

For every k x , k z the solution of (fT4|) yields possible eigen- 
frequencies LOi of the solutions of Maxwell equations that 
should be substituted into the definition of the Casimir 
energy E — J^i^i/^- These solutions should tend to 
zero for y — > ±oo. The summation over the eigenfre- 
quencies is performed by making use of the argument 
principle (fTTj) , which yields the Casimir energy of two 
parallel gratings on a "unit cell" of period d and unit 
length in z direction: 

hed f + °° f +oc f^- 
E = / duj / dk z / dk x 

(2tt) 3 J J_ oq J q 

lndetfj- Ri(iuj)R 2up (iu!)j, (15) 

c is the speed of light in vacuum. This is an exact expres- 
sion valid for two arbitrary periodic dielectric gratings 
separated by a vacuum slit. It can be applied to calcu- 
late the Casimir energy of any parallel periodic gratings 
made of a material described by a dielectric function, 
with surface corrugations of arbitrary geometry. 

Consider the particular case s = 0, depicted in FigfT] 
From the derivation sketched above it follows that 

R 2 -u, P {iu) = K(iu)R 2 (iu)K(iuj), (16) 
where K{iu>) is a diagonal 2(27V + 1) matrix of the form: 

*M = (o g) > < 17 > 

with matrix elements e - L ^" 2 + k l+V c *+ 2j T L )* ( TO 

Y . . . N) on a main diagonal of a matrix G. Note that 
in all Rayleigh expansions the Fourier basis is taken sym- 
metrically around m = 0. When changing the maximum 
value of m from N — 1 to N, each Rayleigh coefficient 
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FIG. 2: Casimir force normalized by its PFA value for two Si 
gratings with a = lOOnm and di = | as a function of d at a 
fixed distance L — 250nm. 

Rpfp(iu)) appearing in the reflection matrices is multi- 
plied by a factor ~ e -^NL/d comm g f rom the matrix 
K(iuj). As a consequence, when 2irNL/d ^ 1 is satisfied, 
the contribution of the coefficients i?Ar p (iw) is suppressed 
exponentially. Therefore for large enough N changing N 
has only a little impact on the final result. 



RECTANGULAR GRATINGS 

We have numerically calculated the exact Casimir force 
for two rectangular gratings at zero temperature in the 
geometry of FigQ] for silicon for different values of d, 
dl = d/2 and a = 100 nm by making use of the for- 
mulas (fT5l [TBI fTTj) and a Drude-Lorentz model for the 
dielectric permittivity of intrinsic silicon 0. We com- 
pare our exact results of the Casimir force for differ- 
ent values of d to the PFA results. Calculated with 
the proximity force approximation, the Casimir force 
between the two gratings is just the geometric sum of 
two contributions corresponding to the Casimir force be- 
tween two plates Fpp at distances L and L — 2a, that 
is Fpfa = \{Fpp{L) + Fpp(L — 2a)). In particular it is 
independent of the corrugation period d. To assess quan- 
titatively the validity of PFA, we plot the dimensionless 
quantity p = |12 |- The ratio is presented on FigUl 

Exact and PFA results differ for silicon by up to 24 per- 
cents for a corrugation period of lOOnm and the PFA vi- 
olation could thus be demonstrated experimentally. We 
recover the PFA result in two limiting cases, for a van- 
ishing corrugation period and for very large corrugation 
periods. In between the exact result for the Casimir force 
is always smaller than the PFA prediction, in contrast to 
calculations for perfect conductors, where the resulting 
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FIG. 3: Casimir force between a Si grating and a gold plate 
as a function of distance. 

force is always larger than the PFA prediction. 

We will now apply our method to the recent experi- 
ment by Chan et al. [l[ , who measured the Casimir force 
gradient between a silicon grating with nanostructured 
trenches and a gold sphere of radius R = 50/xm. The 
force gradient F PS between a sphere of radius R and a 
plate can be expressed via the force Fpp in the plane- 
plane configuration as F' PS = 2nRFpp. This is why we 
show in figure [3] the zero temperature result for the abso- 
lute force values evaluated for a grating with the experi- 
mental parameters a — 980nm, d — 400nm, d\ = 196nm 
placed in front of a gold plate (we used a plasma model 
with a plasma frequency lu p = 9eV for gold and a Drude- 
Lorentz model for intrinsic silicon Q). 

From our calculation we obtain a force Fpp = 
0.51N/m 2 for a plate separation of 150nm. With the ex- 
perimental parameters this leads to a prediction for the 
Casimir force gradient of F' = 160.8,56.4,24.6 pN//xm 
at respectively L — a = 150, 200, 250 nm. The absolute 
values of the force are thus in good agreement with the 
measured values depicted in Fig. 3c of [l|. 

We finally present ratios of our results for the force to 
the predictions of PFA for two different gratings. Figure 
[4] shows p as a function of L — a for two gratings corre- 
sponding to the experiment with a = 980nm, d — 400nm, 
d\ = 196nm (green line) and a = 1070nm, d = lOOOnm, 
d\ = 522nm (blue line) and gives reasonable agreement 
with experimental points and the fit in Fig. 3d of 0]. 

The fact that the perfect conductor model fails might 
be due to the influence of surface plasmons, as the grat- 
ing affects their dispersion relation. Surface plasmons 
contribute essentially and at all distances to the Casimir 
force [HI, 0, El, El, the Casimir force thus has to 
change considerably when structured surfaces are consid- 
ered. These changes are not visible in a perfect conductor 
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FIG. 4: Casimir force normalized by the PFA value between 
a Si grating and a gold plate as a function of distance for two 
different gratings. Solid curves are calculated by making use 
of the least square method from the theoretical points on the 
figure. 



model which ignores the existence of surface plasmons. 
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